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ABSTRACT 

We investigate the evolution of the properties of model populations of ultraluminous X-ray sources (ULXs), 
consisting of a black-hole accretor in a binary with a donor star. We have computed models corresponding 
to three different populations of black-hole binaries, motivated by our previous studies. Two of the models 
invoke stellar-mass (~ IQMq) black-hole binaries, generated with a binary population synthesis code, and 
the third model utilizes intermediate-mass IOOOMq) black-hole accretors (IMBHs). For each of the three 
populations, we computed 30,000 binary evolution sequences using a full Henyey stellar evolution code. A 
scheme for calculating the optical flux from ULXs by including the reprocessed X-ray irradiation by, and the 
intrinsic viscous energy generation in, the accretion disk, as well as the optical flux from the donor star, is 
discussed. We present color-magnitude diagrams (CMDs) as "probability images" for the binaries as well as 
for the donor stars alone. "Probability images" in the plane of orbital period and X-ray luminosity are also 
computed. We show how a population of luminous X-ray sources in a cluster of stars evolves with time. 
The most probable ULX system parameters correspond to high-mass donors (of initial mass > 25 Mq) with 
effective O through late B spectral types and equivalent luminosity classes of IV or V. We also find the most 
probable orbital periods of these systems to lie between 1-10 days. Estimates of the numbers of ULXs in a 
typical galaxy as a function of X-ray luminosity are also presented. From these studies we conclude that if 
the steUar-mass black-hole binaries are allowed to have super-Eddington limited X-ray luminosities: (i) the 
value of the binding energy parameter for the stellar envelope of the progenitor to the black-hole accretor must 
be in the range of 0.01 < A < 0.03 in order not to overproduce the ULXs, and (ii) the stellar-mass black- 
hole models still have a moderately difficult time explaining the observed ULX positions in the CMD. Other 
possible explanations for the apparent overproduction of very luminous X-ray sources in the case of stellar- 
mass black-hole accretors are discussed. Our model CMDs are compared with six ULX counterparts that have 
been discussed in the literature. The observed systems seem more closely related to model systems with very 
high-mass donors in binaries with IMBH accretors. We find that a significant contribution to the optical flux 
from the IMBH systems comes from intrinsic accretion disk radiation whose source is viscous dissipation of 
gravitational potential energy. In effect, the IMBH systems, when operating at their maximum luminosities 
(lO"*' - 10"*^ ergs s~'), are milli-AGN. With regard to the IMBH scenario, while attractive from the aspect of 
binary evolution models, it leaves open the larger question of how the IMBHs form, and how they capture 
massive stellar companions into just the correct orbits. 

Subject headings: (stars:) binaries: general — X-rays: binaries — galaxies: star clusters — accretion, accretion 
disks — black hole physics 



1. INTRODUCTION 

An ultraluminous X-ray source (ULX) is an off-nucleus 
point-like source whose X-ray luminosity, Lj,, exceeds 2 x 
lO''^ ergs s~'. An obvious candidate for such a source is an 
X-ray binary where a compact object accretes from a donor 
star. Under normal circumstances, the accretion rate is con- 
strained by the Eddington hmit of the accretor. The Eddington 
Umit for a neutron star is ~ 10^^ ergs s"^ and that for a IOMq 
black hole (BH) is ~ 10^' ergs s"'. Therefore, if the accretors 
are assumed to be stellar mass BHs, the observed luminosi- 
ties of ULXs exceed the Eddington limit by factors of up to 
~ 100. The question as to how the Eddington limit could be vi- 
olated remains a subject of considerable debate. On the other 
hand, if we assume the accretor to be an intermediate-mass 
black hole of ~ 10^ — W^ Mq, the ULX luminosities can all 
be accounted for without violating the Eddington limit. If the 
accreting stars are indeed IMBHs, the question remains as to 
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how these objects are formed and how they acquire a massive 
stellar companion. 

Several scenarios have been put forth to explain the ap- 
parent super-Eddington luminosities, assuming stellar-mass 
BHs. King et al. (2001) suggested that geometrical beaming 
in the direction of the observer, due to a thick accretion disk 
(Jaroszyriski, Abramowicz, & Paczynski 1980), could lead to 
both a high luminosity and an even higher inferred luminosity. 
Kording et al. (2002) make a similar argument with relativis- 
tic beaming, due to jets, as the explanation. In both of these 
scenarios, the emission is considered to be anisotropic. These 
somewhat contradict the observations of ionization nebulae 
around some ULXs (Pakull & Mirioni 2003) which tend to 
indicate the full implied isotropic luminosity. Another sce- 
nario proposed by Begelman (2002; 2006) suggests that pho- 
ton bubble instabilities in the accretion disk can lead to an 
isotropic luminosity exceeding the Eddington limit by a fac- 
tor of ^10. Socrates & Davis (2006) invoke a hot, optically 
thin corona in conjunction with a geometrically thin, but opti- 
cally thick, accretion disk to explain the observed ULX lumi- 
nosities. Yet other scenarios invoke "shm" accretion disks to 
produce approximately isotropic luminosities that may exceed 
the Eddington limit by up to factors of ~10 (Abramowicz et 
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al. 1988; Ebisawa et al. 2003). Although these scenarios have 
not been shown to work conclusively, they could, in principle, 
account for ULX luminosities up to ~10'*'' ergs s"'. In addi- 
tion to these specific suggestions for bypassing the Eddington 
limit, the continuity and simplicity of the luminosity function 
of luminous X-ray sources from 10^^ — 10^" ergs s~' (a simple 
power law; Grimm et al. 2003) has led some to conclude that 
these represent a single class of systems, i.e., neutron stars and 
stellar-mass black holes^ accreting from a normal donor star. 
However, we note that essentially none of these ideas involv- 
ing stellar-mass BHs would plausibly be able to explain the 
ULXs at the higher end of the observed luminosities, namely 
those with > 10'*" ergs s~'. 

By contrast, the IMBH scenario, as suggested first by Col- 
bert & Mushotzky (1999), accounts for the luminosities of 
all the ULXs because the Eddington limit for an IMBH in 
the mass range ^ lO^-lO'^M© varies between lO'^^'-lO'*^ 
ergs s"'. Other evidence that may suggest the presence of 
an IMBH includes observations of mHz QPOs (Strohmayer 
& Mushotzky 2003) and ionization nebulae (PakuU & Miri- 
oni 2003), and inferences of cool inner accretion disk tem- 
peratures (Miller et al. 2003; Miller et al. 2004; Cropper et 
al. 2004). Despite the several pieces of evidence supporting 
IMBH binaries as the model for ULXs, one of the main prob- 
lems with this scenario Hes in the unknown formation mech- 
anism for such binaries. A few scenarios for IMBH forma- 
tion have been proposed in the literature (see, for example, 
Portegies Zwart et al. 2004; Tutukov & Fedorova 2005) but 
a definitive scenario remains elusive. In particular, the evo- 
lution of supermassive stars with M :» 100 Mq is not well 
understood. 

In three of our earlier papers, we studied the evolution 
of populations of binary systems pertaining to both mod- 
els for ULXs, namely: (i) super-Eddington accretion onto 
LMBHs (Podsiadlowski et al. 2003, Rappaport et al. 2005); 
and, (ii) sub-Eddington accretion onto IMBHs (Madhusud- 
han et al. 2006). In those studies, the main objective was to 
investigate the evolution of the X-ray luminosities of the sys- 
tems with time, and to estimate the formation efficiencies for 
each scenario. In Podsiadlowski et al. (2003) and Rappaport 
et al. (2005), we considered LMBHs in binary systems with 
donor stars in the mass range 2- 17 Mq. Some of the binary 
populations were generated using the binary population syn- 
thesis (BPS) code developed in Podsiadlowski et al. (2003). 
There, we showed that by allowing for a violation of the Ed- 
dington Umit by a factor of ^10-30, such systems would be 
able to explain most of the observed ULXs, except for the 
very most luminous ones. In Madhusudhan et al. (2006), we 
considered donor stars in binary systems with IMBHs with a 
representative mass of 1000 M©. There, we showed that in 
order to have a plausible formation efficiency for ULXs with 
Lx ;J 10"*" ergs s"' the donor stars should be massive (> 8 Mq) 
and the initial orbital separations should be close (< 6-40 
times the radius of the donor star when on the main sequence). 

A number of fairly secure optical identifications of ULX 
counterparts have now been reported in the literature. Liu et 
al. (2002) identified a unique counterpart to the ULX NGC 
3031 XI 1 (in M81) and found it to be consistent with an 08 
main sequence (MS) donor star. Kaaret et al. (2004) reported 
an optical counterpart to the ULX in Holmberg II to be con- 

In this paper, we refer to a stellar-mass black holes as a "low-mass black 
hole" or "LMBH", rather than "SMBH" which might be confused with "su- 
permassive black hole". 



TABLE 1 

Magnitudes and Colors of ULX counterparts 

(ADAPTED FROM COPPERWHEAT ET AL. 2007) 



System 


Mb 


Mv 


B-V 


NGC 4559 X-7 (CI) 


-7.22 ±0.19 


-7.03 ±0.16 


-0.19 ±0.25 


M81 X-6 


-4.28 ±0.04 


-4.18 ±0.03 


-0.10 ±0.05 


MIDI ULX-1 


-6.19±0.15 


-5.92 ±0.12 


-0.27 ±0.19 


NGC 5408 ULX 


-6.40 ±0.20 


-6.40 ±0.20 


-0.00 ±0.28 


Holmberg H ULX 


-6.03 ±0.19 


-5.78 ±0.11 


-0.25 ±0.22 


NGC 1313 X-2 (CI) 


-4.70 ±0.18 


^.50±0.18 


-0.20 ±0.25 



sistent with a donor star with spectral type between 04 V and 
B3 lb. Kuntz et al. (2005) have studied the optical counterpart 
of M 101 ULX-1 and found the colors to be consistent with 
those for a mid-B supergiant. Soria et al. (2005) studied can- 
didate counterparts to the ULX in NGC 4559 and suggested 
the donor to be either a blue or a red supergiant of high mass 
(~1O-2OM0). Mucciarelli et al. (2005; 2007) have found 
candidate optical counterparts to the ULX NGC 1313 X-2 that 
are either BO-09 main-sequence stars or G supergiants. 

The best candidate optical counterparts to the ULXs dis- 
cussed above have been conveniently selected and summa- 
rized by Copperwheat et al. (2007). In particular, they tabu- 
late the photometric values for the various systems, calculat- 
ing reddening corrections and absolute magnitudes, wherever 
necessary. As can be seen from the above list, the observa- 
tions seem to indicate that the spectral classes of the most 
promising candidates generally range between O and B spec- 
tral type, and the luminosity classes are either V or lb. In the 
present study, we consider optical counterparts of six ULX 
systems for which photometric data in the B and V bands are 
available in the Uterature. The photometric data for these sys- 
tems are given in Table 1 (as adapted from Copperwheat et 
al. 2007). 

Theoretical studies of ULXs in the optical regime have 
heretofore concentrated on constraining the nature of the 
donor star, the mass of the accretor, and the orbital period. 
All the models follow the standard paradigm of ULXs be- 
ing X-ray binaries with active accretion through Roche -lobe 
overflow from a donor star onto a BH accretor. Rappaport 
et al. (2005) discussed preliminary theoretical models involv- 
ing disk irradiation in ULXs. They presented sample evo- 
lution tracks, on a color-magnitude diagram, for four ULX 
models, and discussed the optical appearance of ULX BH bi- 
naries. Pooley & Rappaport (2005) suggested detection of 
X-ray and/or optical eclipses as a means to constrain the mass 
of the accretor. In recent studies, Copperwheat et al. (2005; 
2007) have also constructed irradiation models to describe 
the optical emission from ULXs, and used the models to con- 
strain the properties of several systems observed in the opti- 
cal. The parameters being considered were the mass, radius 
and age of the donor, and the BH mass in some cases. From 
the fits to the various observations, they find the counterparts 
to be consistent with being main-sequence stars or evolved 
giants/supergiants with spectral types O, B, or A. 

In the present paper, we report a detailed population study 
spanning a large region of parameter space of ULX properties 
in an effort to better explore and help constrain ULX models. 
In particular, we investigate models of the optical properties 
of ULXs. We choose three sets of representative populations 
from our previous studies, and explore their optical proper- 
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ties and other observables. For the models with LMBHs, we 
choose two sets of populations obtained from the BPS code 
corresponding to two different prescriptions for the critical 
mass ratio in the primordial progenitor binary that determines 
whether a common envelope phase can occur. These two 
models yield very different distributions of secondary (i.e., 
ULX donor) star masses. For models with IMBHs, we choose 
model C from Madhusudhan et al. (2006). This model yielded 
the highest formation efficiency for ULXs in the IMBH sce- 
nario. For, each of the three populations, we follow 30,000 
binary evolution calculations. The models are summarized in 
Table 2. 

The X-ray luminosities are calculated using a standard for- 
mulation of Roche-lobe mass transfer in X-ray binaries. The 
optical flux from the system is determined as the sum-total of 
the optical flux from the donor star, X-ray irradiation of the 
accretion disk, and intrinsic energy generation in the disk. We 
present detailed color-magnitude diagrams (CMDs) for all the 
models. These include CMDs for the binary system, i.e., the 
sum of the optical flux of the donor and that due to radiation 
from the disk, as well as CMDs for the donor star alone. We 
also study in detail the evolution of the X-ray luminosity with 
the age of the system, fev> and the evolution of the X-ray lu- 
minosity with orbital period, Porb- 

2. METHODS 

2.1. Binary Population Synthesis 

In our previous studies we explored several models of 
ULX populations for both the low- (hereafter "LMBH") and 
intermediate-mass (hereafter "IMBH") black-hole scenarios 
(Rappaport et al.2005; Madhusudhan et al.2006). In the 
present work, we chose two LMBH- and one IMBH-binary 
population that were motivated by these previous studies. The 
parameters for the three models are suimnarized in Table 2, 
and the distributions of the initial system masses and orbital 
periods are shown in Figure 1. Models La and Lb are LMBH 
populations, whereas Model Ic represents an IMBH popula- 
tion. 

To generate the LMBH binary populations, we used the bi- 
nary population synthesis (BPS) code developed in Podsiad- 
lowski et al. (2003). We briefly review here the formulation 
from Rappaport et al. (2005). We started with a very large set 
of massive primordial binaries and generated a much smaller 
subset of these that evolved to contain a black hole and a rela- 
tively unevolved companion star. The product was a set of 
"incipient" black-hole X-ray binaries with a particular dis- 
tribution of orbital periods, foi-b, donor masses, Mjon, and 
black-hole masses, Mbh, for each of a number of different 
sets of input assumptions (see, e.g.. Fig. 2 of Podsiadlowski 
et al. 2003). For this part of the calculation, we employed 
various "prescriptions", based on single-star evolution mod- 
els for the primary, simple orbital dynamics associated with 
wind mass loss and transfer, assumptions about the magni- 
tude of the wind mass loss from the primary as well as from 
the core of the primary after the common envelope, and na- 
tal kicks during the core collapse and formation of the black 
hole. 

The decision in the binary population synthesis code of 
whether or not a common-envelope phase occurs, when the 
primary first overflows its Roche lobe and starts to transfer 
mass onto the secondary, is based on the evolutionary state of 
the primary when mass transfer commences. We first define 
^TAMS and Ryig as the radii of the primary when it reaches 
the terminal main sequence and the end of the Hertzsprung 



TABLE 2 
ULX Population Models 



Model^ 


Mbh'' 






A= 


n ■ f 
ycnt 


La 


6-15 


5-18 


1-18 


0.1 


1.2 or 2.0 


Lb 


6-15 


5-27 


1-12 


0.1 


1.2-^^2.0 


Ic 


1000 


5-50 


1-10 







"The distributions of masses and periods for the different models are shown 
in Figure 1 

''Approximate range of black hole masses, in units of Mq 
'' Approximate range of donor masses, in units of Mq 
''Approximate range of initial orbital periods, in units of days 
^Dimensionless inverse binding energy of the envelope of the black-hole 
progenitor star 

*For details of the prescription see eqns. (1) - (3) in the text. 



gap, respectively. Given the initial orbital separation and the 
mass ratio of the primordial binary, q = Mprim/Msec, we can 
compute the Roche -lobe radius, Ri, of the primary. In or- 
der to form a relatively close black-hole binary of the type 
we are considering in this work (i.e., /orb ^ 20 days), we re- 
quire that the initial mass-transfer rate from the primary to 
the secondary be so large that it leads to a common-envelope 
and spiral-in phase. This depends on the initial mass ratio and 
evolutionary state of the secondary and is somewhat uncertain 
(see Pols 1994; WeUstein, Langer & Braun 2001). For one of 
our two models involving stellar-mass black hole accretors, 
we utilize the follow critical mass ratios, ^crit> depending on 
the relation among Rl, /?tams, and Rhc'- 

^crit = 2 for Rtams <Rl< Rhg (1) 
gcrit=l-2 for Rl> Rhg (2) 

For the second of our two models involving stellar-mass 
black-hole accretors we utilize an interpolated version of the 
above prescription: 



^crit = 2.0-0.8 



l0g/?L-l0g/?TAMS 



log/?HG-log/?TAMS 



(3) 



for /?TAMS < Rl < Rhc- This somewhat ad hoc prescription 
has the net effect of producing a much greater fraction of in- 
cipient ULX sources with high-mass donor stars. 

Simple energetic arguments were used to yield the final-to- 
initial orbital separation during the common-envelope phase 
wherein the envelope of the primary is ejected. Here we uti- 
lized a parameter A, which is the inverse of the binding energy 
of the primary envelope at the onset of the common-envelope 
phase in units of GMpnjaMg/Rprim, where M^^, Mg and Rp^m 
are the total mass, envelope mass, and radius of the primary, 
respectively. This parameter strongly affects the final orbital 
separation after the common-envelope phase, where smaller 
values of A correspond to more tightly bound envelopes, and 
hence more compact post-conmion-envelope orbits. 

With the above parameterization for the ejection of the 
common envelope we find the foUowing expression for initial- 
final orbital separation: 



CE 



M, 



pnm 



2Me 



(4) 



(e.g., Webbink 1985; Dewi & Tauris 2000; Pfahl et al. 2003), 
where the subscripts "prim", "c", and "e" stand for the pro- 
genitor of the black hole, its core, and its envelope, respec- 
tively, and "sec" is for the secondary, i.e., the progenitor of 
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the "donor star" in the black-hole binary. The quantity is 
the Roche-lobe radius of the black-hole progenitor in units 
of a,, 8cE is the fraction of the gravitational binding energy 
between the secondary and the core of the black-hole pro- 
genitor that is used to eject the common envelope, and A 
is defined above. For typically adopted parameter values, 
A ~ 0.01- 1 (e.g., Dewi & Tauris 2000; Podsiadlowski et al. 
2003), ScE — 1, and ~ 0.45-0.6 (for an assumed mass ratio 
between the black-hole progenitor and the companion in the 
range of~2:1^15:l), the second term within the paren- 
theses in eq. (4) dominates over the first. In this case, we find 
the following simpUfied expression for a//fl,: 

where the leading factor is r/,/2 ~ 1/4, while the factor in 
parentheses involving the black-hole progenitor is ^ 0.020 ± 
0.002 Mq for virtually all of the progenitors we consider. 
This explains why the large majority of the incipient black- 
hole binaries (with low- to intermediate-mass donors) found 
by Podsiadlowski et al. (2003) resulted from an initially very 
wide orbit (Poib ~ years - when the primary attains radii of 
^1000-2300 Rq) preceding the common-envelope phase in 
order to avoid a merger between the secondary and the core 
of the primary. 

In the present work, for Models La and Lb (see Table 2 
and Fig. 1), we chose two populations that were obtained us- 
ing the BPS code described above. These two LMBH pop- 
ulations both utilize an inverse envelope binding energy pa- 
rameter A = 0.1, and two different prescriptions for the criti- 
cal mass ratio required to produce a common envelope phase 
(see eqns. [l]-[3]). These result in quite different produc- 
tion rates for incipient ULXs with high-mass donor stars (i.e.. 

For the IMBH case. Model Ic (see Table 2), we chose the 
same population as Model C in Madhusudhan et al. (2006). 
Due to our lack of a firm understanding of how IMBHs form 
in star clusters, and how they capture companion stars, we 
adopted a very simple prescription for the population of in- 
cipient IMBH binaries (but, see Blecha et al. 2006). For each 
binary, we chose the initial donor mass to lie uniformly in the 
range of 5-50 Mq using Monte Carlo methods. And, we 
chose the initial orbital separation uniformly over the range 
of 1 - 3 times the separation required for a ZAMS star of the 
corresponding mass to fill its Roche lobe. This translates to a 
range of initial orbital separations between 40-200 Rq, and 
Porb in the range of ~1 - 10 days. We found in our previ- 
ous work that this population is the most favorable one we 
explored for producing ULXs in the IMBH scenario, i.e., the 
donor stars should be massive, i.e., > 8 Mq and the initial 
orbital separations, after circularization, should be close, ^ 
6-40 times the radius of the donor star when on the ZAMS. 
This range of close orbital separations represents the case of 
direct tidal capture and circularization of a single field star by 
the IMBH (see, e.g., Hopman et al. 2004; Pfahl 2005). 

For, each of the three modeled ULX populations, we carried 
out 30,000 binary evolution calculations. All models were run 
using 60 nodes of the elix3 Beowulf cluster located at the 
University of Sherbrooke, Quebec. The run time for each of 
the three models was ^^40 hours. 

2.2. Stellar Evolution 

The stellar evolution of the donor stars, including mass 
loss, was followed with EZ which is a stripped down, rewrit- 
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Fig. 1. — Distributions of initial donor masses, black hole masses and 
initial orbital periods of the LMBH systems. The red curves correspond to 
Model La and the green curves coirespond to Model Lb. In the top panel, 
the solid curves represent donor masses and the dashed curves represent BH 
masses. For Model Ic, Mbh = 1000 Mq and we choose M^an uniformly 
between 5-50 Mq, and the initial orbital separation uniformly between 1-3 
times the separation required for a ZAMS star of the corresponding donor 
mass to fill its Roche lobe. 

ten version of a subset of the stellar evolution code devel- 
oped by P. P. Eggleton (Paxton 2004). The physics of the 
program is unchanged from Eggleton's (essentially as de- 
scribed in Pols et al. 1995), but the structure of the code 
has been modified to facilitate experiments involving pro- 
grammed control of parameters. There are zero-age main- 
sequence (ZAMS) starting models for a variety of metallici- 
ties (from Z = lO""* to Z = 0.03) and masses (from 0.1 to 100 
Mq), with arbitrary starting masses created by interpolation. 
A user-provided procedure is called between steps of the evo- 
lution to inspect the current state, to make changes in param- 
eters, and to decide when and what to record to log files. The 
source code and data for E Z can be downloaded from the web 
at <http://theory.kitp. ucsb.edu/^paxton>. For all models in 
this particular study, the number of stellar mesh-points was 
fixed at 200 in the interest of minimizing computation time. 

2.3. Binary Evolution Calculations 

The binary evolution was governed by a sequence of cal- 
culations involving the mass-transfer rate, the corresponding 
change in orbital separation, and the subsequent monitoring 
of mass loss. The mass transfer considered was due solely to 
Roche-lobe overflow. When the donor star fills its Roche lobe, 
the excess matter flows through the inner Lagrange point onto 
the accretor. Assuming spherical geometry, the mass-transfer 
rate was calculated using M ~ lirRMpv, where R is the radius 
of the donor, H is the density scale height of the donor atmo- 
sphere, p is the density of the atmosphere at the Roche lobe, 
and V is the thermal velocity at the photosphere. Under the 
approximation of an isothermal, constant-gravity atmosphere, 
the scale height is given by // ~ A:r/ pg, where p is the mean 
molecular weight, and the density profile of the atmosphere is 
exponential. The Roche-lobe radius of the donor was taken to 
be/;RL = 0.49fl^2/3[o.6^2/3_^ljj(^_^^i/3)]-i (gggigton 1983), 

where, a is the separation of the binary, q = MionlM^a, Mdon is 
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the mass of the donor, and Mbh is the mass of the black hole. 
This procedure, while making use of some approximations, is 
self-adjusting to yield the correct mass-transfer rates. 

The Eddington limited mass-transfer rate onto the accre- 
tor is given by MEdd = ^i^GMacc/'H^c, where G is the gravi- 
tational constant, k is the radiative opacity, and rj is the ef- 
ficiency of the black-hole accretor in converting rest mass to 
radiant energy. The opacity is assumed to be predominantly 
due to electron scattering and is given by k = 0.2(1 -l-Z) cm^ 
g~\ where, X is the hydrogen mass fraction. For a spinning 
black-hole, we take the accretion efficiency to be given by 
?7 = 1 - a/1-(M/3Mo), where Mq is the initial mass of the 
black-hole (Bardeen 1970). This assumes an initially non- 
spinning black-hole. 

If we consider all the matter leaving the donor to be retained 
by the accretor, independent of the Eddington limit, then the 
resulting luminosity is referred to as the "potential luminos- 
ity" and is given by, Lpot = rjMc^. The actual isotropic X-ray 
luminosity, on the other hand, is given by = PrjMc^, where 
(3 restricts the mass-transfer rate to the Eddington limit. For 

M < MEdd, = 1. Otherwise /3 =MEdd/M- When M > Mndd. it 
is assumed that all the matter passes through the accretion disk 
until it reaches the inner edge of the disk, at which point the 
excess matter (above the Eddington hmit) is ejected out of the 
system in the form of a jet. The mass and orbital angular mo- 
mentum lost from the system in such ejection is incorporated 
when calculating the orbital separation of the system. We em- 
phasize that we use the Eddington limited mass-transfer rate 
while calculating the orbital parameters during the evolution 
of the binary. However, along with the various orbital param- 
eters, we also record Lpot at each step. And, in all the results 
presented in this paper, we typically refer to Lpot instead of 
the X-ray luminosity, Li, of the system in order to examine 
the potentialities of violating the Eddington limit. 

2.4. Disk Irradiation 

The optical flux from a ULX is comprised of contributions 
from the donor star and from the reprocessing of X-ray pho- 
tons by the accretion disk (as well as intrinsic energy gener- 
ation within the disk). The irradiation of the accretion disk 
depends upon L^, the geometrical properties of the disk, and 
whether or not the central accreting star has a hard surface. 
For purposes of the current calculation, we assume the disk 
to be geometrically thin. We find the effective temperature of 
the disk to be: 



r(r)~ 



1/4 



e'Vio/^(l-a) + 3x- 



1/4 



(6) 



Here, x = r/r^i„, where r is the radial distance, r^i^ is the in- 
ner disk radius, taken to be 6GMbh/c^, ct is the X-ray albedo 
of the disk, which we take to be 0.7. The half thickness of 
the disk is given by h{r) = ir'^l''f^ll, and ^' = C(''min/''max)"''^, 
where ^ is a constant equal to /j('"max)/''max and rmax is the 
outer radius of the disk. We have taken rmax to be O.lrLa 
where r^ is the dimensionless Roche-lobe radius of the black- 
hole accretor. We somewhat arbitrarily adopt a value for ^ 
of 0.1, corresponding to a full angular thickness of the disk 
of ~ 12°. The second term in the square brackets in eq. (6) 
comes from the Shakura-Sunyaev (1973) solutions for the 
temperature profile of a thin disk around a BH accretor. We 
have taken the inner edge of the accretion disk to be at the 
iimermost stable circular orbit around a non-rotating BH. We 



have neglected the factor ( 1 - ■\/ r^in/r) since its contribution 
to the optical flux can be assumed to be negligible. The first 
term in the square brackets in eq. (6) is adapted from the ex- 
pression for the irradiation temperature of an accretion disk 
around a BH derived in King, Kolb & Szuszkiewicz (1997). 
This is a modified version of the corresponding term in eq. (8) 
of Rappaport et al. (2005) where we used an expression ap- 
propriate for radiation emanating from a centrally located, 
isotropicaUy emitting, hard surface. The present formulation 
takes into account the fact that the X-ray radiation comes from 
the innermost part of the accretion disk surface and is emitted 
with a Lambertian angular distribution. Since the irradiated 
portion of the outer disk lies at a shallow angle with respect 
to the surface of the inner disk, the small value of the cosine 
factor significantly reduces the X-ray irradiation. 

As an important caveat to our disk irradiation calculations, 
we note that the assumption of a thin disk (emitting in a Lam- 
bertian manner) likely breaks down at accretion rates slightly 
or greatly exceeding Eddington. In this case, the irradiation 
of the outer disk might become more or less efficient. How- 
ever, such a calculation is considerably beyond the scope of 
the present work. 

The optical flux from the disk is determined by using the 
temperature profile described in eq. (6), and assumes the local 
spectrum in each annulus of the disk is thermal. The radiant 
flux from each annular ring on the disk surface is determined 
from the specific intensity of the blackbody radiation and the 
temperature at that radius. The total flux due to the disk at a 
particular wavelength is calculated by integrating the annular 
flux over the radial extent of the disk, from r^^^ to Tmax- The 
flux from the donor star is determined from its surface tem- 
perature, Te, as obtained from the stellar evolution code. For 
the wavelengths under consideration, a thermal blackbody is a 
reasonable approximation in order to determine the flux from 
the star (with ~10% accuracy). This corresponds to an error 
in the apparent V magnitude of ~ ± 0.1, which is within the 
error bars of most observations, and negligible compared to 
the absolute values of V and My- The corresponding error in 
B-V is < 0.2. 

The total flux of the system consists of contributions from 
the disk and the donor star. The B and V magnitudes were 
calculated from the total fluxes using B = -2.5 logiFs / Fbq) 
and V = -2.5 log(/v/-PVo), where Lg, Fbq, Fy and /vo are the 
fluxes and reference fluxes at A = 4380 A and A = 5450 A, 
respectively {Fbo = 6.61 x 10~' ergs cm~^ s"' A"' and FyQ = 
3.64 X 10-' ergs cm^^ s"! A^i). 

After we compute B-V for the combined radiation from the 
donor star and the accretion disk^, we apply a small correction 
(< 0.2 magnitudes) as a function of B-V to take into account 
the treatment of the radiation as blackbody emission. These 
corrections were based on the differences between the B-V 
colors computed from the simple algorithm given above and 
the colors of main sequence and supergiant stars tabulated by 
Johnson (1966). 

2.5. Intrinsic Disk Emission 

We have found from our models (see §3) that under cer- 
tain circumstances the intrinsic disk emission (due to the vis- 
cous release of gravitational potential energy) which is rep- 
resented by the second term in square brackets in eq. (6), 
can dominate the optical emission from the disk irradiation 

' All of our binary systems are assumed to be viewed at an average orbital 
inclination angle of 60° . 
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Fig. 2. — Illustrative temperature profiles of the irradiated accretion disk 
for an IMBH model (blue curves) and LMBH model (red curves). The solid 
curves are for the case of intrinsic (viscous) heating of the accretion disk 
alone, while the dashed curves are for the iiradiation of the disk without any 
contributions from viscous heating. These correspond to the 2nd and the 
1st terms, respectively, in the square brackets of eq. (6). The binary system 
parameters used for this illustrative example are given in the text. 

(first term in square brackets). To illustrate this effect, we 
show in Fig. 2 the temperature profiles for an illustrative set 
of binary parameters: Lpot = lO'*" ergs s"', Porb = 30 days, and 
Mdon = 10 Mq. For an IMBH model (Mbh = 1000 M©) the 
minimum and maximum disk radii are r^in = 9 x 10** cm and 
'"max = 1.5 X 10'^ cm, respectively, while for an LMBH model 
(with Mbh = 10 Mq), these radii are rmin = 9 x 10^ cm, and 

'"max = 2 X 10'^ cm. 

The red and blue curves in Fig. 2 are for the LMBH and 
IMBH models, respectively. The dashed curves are the tem- 
perature profiles for the case where irradiation is the only 
source of energy input to the disk; conversely, the solid curves 
are for the intrinsic (viscous) release of gravitational energy 
alone - without X-ray irradiation. For the LMBH model, we 
see that the intrinsic energy release dominates the irradiation 
contribution to the optical radiation only for radial distances in 
the disk of < 1 Rq. Since the optical light from the donor star 
itself comes from a typically much larger area, the intrinsic 
disk emission is generally not competitive with that from disk 
irradiation. By contrast the crossing point for the two tem- 
perature profiles corresponding to the IMBH model occurs at 
a radial distance of ~55 Rq (and T(r) remains > 10"* K out 
to > 20 Rq). These radii are larger than the size of the donor 
star while it is near the main sequence or subgiant branch (i.e., 
while it is hot). Thus, it is quite possible that the largest con- 
tribution to the optical light from the IMBH systems comes 
from the intrinsic radiation emitted by the accretion disk it- 
self. We will see the effect of this in §3 (see especially Fig. 6). 

2.6. Irradiation of the Donor Star 

We do not consider the effects of X-ray irradiation of the 
donor star. For our IMBH models the half angle subtended by 
the donor star ranges between ^^5° and 10°. Since we have 
assumed an accretion disk with only a modest half angle of 6°, 
most or all of the donor star can be expected to be shielded by 
the accretion disk. For the LMBH models the donor stars, at 



least initially, subtend larger half angles of ^18° to 24°. The 
fractional solid angle that the donors subtend ranges between 
~0.025 and 0.036 while they are on the main sequence (and 
even less when they are on the giant branch and have trans- 
ferred a significant amount of their mass). After taking into 
account shadowing by the accretion disk, the albedo of the 
donor star surface, and the Lambertian angular dependence 
of the disk irradiator, we find that no more than 0.0018 
0.0026 of the X-ray luminosity is reprocessed on the face of 
the donor star. This is less than comes from the irradiated 
accretion disk and, furthermore, the larger distance between 
the X-ray source and the donor compared to that of the irradi- 
ated disk renders the flux even smaller yet. Therefore, for this 
work, we have neglected the irradiation of the donor stars. 

3. RESULTS AND DISCUSSION 

We present models of ultraluminous X-ray sources includ- 
ing calculations of the optical flux from the system. The op- 
tical flux includes contributions from the donor star as well 
those due to radiation from the accretion disk. As mentioned 
previously, we have chosen three representative populations 
of binary systems for our study. Two of these populations 
have LMBHs and one population has IMBHs as accretors. At 
the beginning of each evolution, the donor star is on the zero 
age main sequence (ZAMS), and the initial system parame- 
ters are set by the binary population synthesis algorithms de- 
scribed above. 

3.1. Evolution of a Single Binary 

Figure 3 shows sample evolutionary tracks for three indi- 
vidual ULX binaries. The three binaries are illustrative of the 
three models listed in Table 2. Starting from the upper left 
panel, panel (a) shows the evolution of Lpot with the age of 
the system (fev). Panel (b) presents the evolution of Lpot with 
the orbital period, Porb- Panel (c) shows the track in the color- 
magnitude diagram (CMD) of the ULX binary, taking into 
account the optical flux from the donor star as well as the flux 
due to radiation from the disk. And, panel (d) represents the 
CMD of the donor star alone. In calculating the optical flux 
from the disk, Lpot was used to irradiate the disk, which allows 
for possible violation of the Eddington limit. The evolution 
of the various parameters with time is denoted by the arrows. 
The blue curves represent the evolution of an IMBH binary 
with an initial donor mass, Mjon = 30 Mq, a black-hole mass, 
Mbh = 1000 Mq, and an initial orbital period, forb = 2.1 days; 
taken from Model Ic (see Table 2). The green curves represent 
Model Lb with Mjon = BMq, Mbh = 12Mq, and Porb = 1-9 
days. And, the red curves are for Model La with Mdon = 1 Mq, 
Mbh = IIMq, and Porb =1.6 days. In order to illustrate the 
direction of evolution of the tracks, on each panel arrows are 
marked on the track that is most clearly separated from the 
others. 

Let us consider the Lpot -fev tracks in panel (a) for illustra- 
tion, and follow the blue curve. The evolution starts with the 
donor on the ZAMS. As the donor evolves through the main 
sequence and somewhat beyond, the radius increases slightly 
and the donor fills its Roche lobe. Mass transfer then takes 
place onto the accretor through the inner Lagrange point. The 
Lpot - tew curve depicts a modest increase in Lpot as the donor 
evolves through the main sequence. For high-mass stars, as 
are considered here, there is an overall contraction phase be- 
fore the star ascends the giant branch. During the contraction 
phase, even a small decrease in radius leads to a large reduc- 
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Fig. 3. — Evolutionary tracks of illustrative binary evolution calculations. Panel (a): evolution of X-ray luminosity with age of the system. Panel (b): Orbital 
period versus X-ray luminosity. Panel (c): Color-magnitude diagram for the binary system. Panel (d): Color-magnitude diagram of the donor star alone. The 
blue curves show the evolution of an IMBH binary, typical of Model Ic. with an initial donor mass (Mjo,,) of 30 Mq, a black hole mass (Mbh) of 1000 Mq, and 
an initial orbital period (Porb). of 2.1 days. The green curves represent Model Lb with Mjon = 13Mq, Mbh = 12M0, and Po,b = 1.9 days. And, the red curves 
represent Model La with Mdon = 7 Mq , Mbh = 1 1 M© , and Porb = 1.6 days. 
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Fig. 4. — Population diagrams for Model La, LMBH with A = 0.1 and t^crit for the primordial binaries taking a value of 1.2 or 2.0 (see text, eqs. (1) and (2), 
and Table 2 for details). Panel (a): Evolution of potential X-ray luminosity with age of the system. Panel (b): Orbital period versus potential X-ray luminosity. 
Panel (c): Color-magnitude diagram of the binary systems. Panel (d): Color-magnitude diagram for donor stars alone. For both panels (c) and (d), contributions 
to the CMD image are made only during times when Lpat ^ 2 X 10^' ergs s"' . Each diagram is an image of 700 X 700 pixels, and represents the 30,000 X-ray 
binary evolution calculations that we computed. For each diagram, the parameter values from all the tracks were registered in each pixel that was traversed. The 
intensity in panel (a) represents the number of systems in the pixel, whereas the intensity in each of the other panels represents the accumulated evolution time 
spent in the pixel. The colors are scaled according to the 1/4 root of the corresponding intensities, with purple being of highest intensity and red being the lowest; 
the actual ratio of values between purple and red is '~200: 1. In panel (c) the irradiation is taken to be from the full Lpm. 



tion in the mass-transfer rate, M, and hence the observed dip 
in Lpot. After the contraction phase, the star ascends the giant 
branch expanding through the Hertzsprung gap on a thermal 
timescale. This increase in radius leads to a spike in Lpof The 
effects of the various phases of evolution of the donor star 
also manifest themselves through the other properties of the 
system. For example, the development of Porb through the var- 
ious evolutionary phases is evident from panel (b). One can 
see the sharp increase in the orbital period as the star expands 



on the giant branch. The relatively smaller variations in Lpot 
as the period changes indicate the near constant, even if rapid, 
rate of growth of the donor star. Panels (c) and (d) depict the 
evolution of the optical properties of the system. As is appar- 
ent from panel (d), the increase in the absolute magnitude of 
the donor star as it ascends the giant branch takes place along 
with a corresponding decrease in its effective temperature. In 
the context of a binary system, however, some of the X-ray 
flux is reprocessed in the accretion disk and the tracks reflect 
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Fig. 5. — Population diagrams for Model Lb, LMBH with A = 0. 1 and for the primordial binaries taking a value over the range of 1.2—2.0 (see text, eq. (3), 
and Table 2 for details). Panel (a): Evolution of potential X-ray luminosity with age of the system. Panel (b): Orbital period vs. potential X-ray luminosity. Panel 
(c): Color-magnitude diagram of the binary systems. Panel (d): Color-magnitude diagram for donor stars alone. All other descriptors are the same as for Fig. 4. 



the effects of mass transfer, as seen from panel (c). 

3.2. Population Diagrams 

Figures 4, 5 and 6 show color images of 30,000 evolution- 
ary tracks each, corresponding to Models La, Lb and Ic, re- 
spectively (see Table 2). The panel arrangement and parame- 
ters shown are the same for all three figures, and match those 
of Fig. 3. Each panel contains a 700 x 700 image matrix of 
the corresponding parameters. For panel (a), the matrix cov- 
ers 7 decades in Lpot and 3.5 decades in evolution time, in 
equally spaced logarithmic intervals. For each step of an evo- 
lutionary track, the position of the Lpot ~ fev P^ir is located in 
the matrix and a value of "1" is added to the corresponding 
matrix element. The resulting matrix, after recording all the 



evolution steps from all the 30,000 tracks, is displayed as a 
color image, with 1/4-root scaling in intensity to enhance the 
dynamic range. For panel (b), the matrix covers 4 decades in 
Porb ™d 7 decades in Lp^t, in equally spaced logarithmic inter- 
vals. Each time the step of an evolution track lands in a matrix 
element, the evolution time-step is added to the corresponding 
matrix element, as opposed to adding a value of ' 1 ' as was the 
case for panel (a). For panels (c) and (d), the matrices cover 
1 1 units in the absolute visual magnitude. My, and 2 units in 
color index, in equally spaced linear intervals. Similar to the 
Porb versus Lpot image in panel (b), a value equal to the evo- 
lution time step is added to the CMD matrix each time the 
step of an evolution track lands in a particular matrix element. 




Fig. 6. — Population diagrams for IMBH Model Ic: Panel (a): Evolution of potential X-ray luminosity with age of the system. Panel (b): Orbital period versus 
potential X-ray luminosity. Panel (c): Color-magnitude diagram of the binary systems. Panel (d): Color-magnitude diagram for donor stars alone. All other 
descriptors are the same as for Fig. 4. 



However, it is extremely important to note that contributions 
to the CMD are made only during times when ^ 2 x 10-''^ 
ergs s~' i.e., when the source would be a potential ULX. 

As in panel (a), the matrices in the remaining three pan- 
els are also displayed as color images with 1/4-root scaling 
in intensity (i.e., relative probability) to enhance the dynamic 
range. The colors reflect this scaling, violet being of highest 
intensity and red being the lowest. 

Thus, the intensity in panel (a) is a measure of the number 
of systems with a particular Lpot found at a particular time in 
the evolution of the cluster. These numbers allow one to cal- 
culate the number of active ULX systems in the population 
at any given epoch and to make estimates of the numbers of 



ULXs of a certain Lpot in typical galaxies. For all the other im- 
ages, the intensity at a point is a measure of the total duration 
(time) spent by all the systems in that interval of parameter 
space. Since this intensity incorporates both the numbers of 
systems and the amount of time each system spends in that in- 
terval of parameter space, it provides a probability map of the 
corresponding parameter space. For purposes of discussion 
later, we shall refer to violet and dark blue regions in these 
diagrams as regions of high probability and to regions in red 
as those of low probability. 

3.3. Evolutionary Phases of the Donor Stars 

Comparing panel (c) from Fig. 3, to that of Figs. 4, 5, and 
6, we see that the regions of high probability correspond to 
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the early phases of the donor star's evolution, i.e., on the main 
sequence or the sub-giant branch phase. This is apparent be- 
cause the donor star spends a predominant amount of its Ufe- 
time on the main sequence, and evolves rather quickly on the 
giant branch. Consequently, the regions of low probability on 
the diagram are the regions corresponding to the giant-branch 
phase of the evolutionary tracks or other regions not readily 
accessed during the binary evolution. However, this descrip- 
tion of the phases of evolution with respect to regions on the 
CMD is only approximate. The exact structure of the CMD 
is a result of a complex interplay between the evolution of the 
donor star, its mass loss, and the radiation from the disk at 
a particular evolutionary stage. From panel (c) of Figs. 4, 5, 
and 6, we find that, with a high probability, the donor stars 
belong to the equivalent luminosity classes for single stars of 
IV and V, with aB-V color range of about -0.35 to -0.10 
(corresponding to an effective spectral class of O through late 
B). 

In many of the panels in Figs. 4, 5, and 6 one can discern 
some of the individual tracks, roughly parallel to each other 
in the initial stages of the evolution. Tracks higher on the 
diagram (i.e., lower My) correspond to higher-mass donors 
since, at any evolution phase, more massive stars are brighter. 
Higher-mass donor stars become progressively rarer and also 
spend lesser amounts of time on the main sequence; hence the 
probability fades toward lower My. 

3.4. Relevance to the Observations 

The data from the optical observations of six ULXs dis- 
cussed earlier (see Table 1) are plotted along with the model 
CMDs. For models La (Fig. 4) and Lb (Fig. 5), we see that two 
of the six data points fall on high-probability regions of the 
CMD (panel c), two on regions of moderate probability, and 
two on low-probability regions. However, for the IMBH sce- 
nario all six data points fall on or close to the high-probability 
regions of the CMD in Fig. 6. Thus, in the framework of the 
models presented here, the IMBH scenario appears to be the 
more favorable one. 

We also see from Fig. 6 that all the observational points fall 
in regions corresponding to the very high donor masses. This 
indicates how the IMBH models fit the data better than the 
LMBH models. The brighter magnitudes of the observed sys- 
tems require massive donor stars. And, massive donor stars 
need more massive accretors in order for mass transfer to oc- 
cur stably via Roche-lobe overflow. In the LMBH scenario, 
the donor masses are limited by the fact that they cannot be 
more than about twice the mass of the BH accretor in order 
for stable mass transfer to take place. In practice, this lim- 
its most of the donor masses to be < 20 Mq. In the IMBH 
case, on the other hand, the BH mass has been set at lOOOM© 
and the maximum donor mass we consider is 50 Mq. Quite 
clearly, it is the higher- mass donors (M > 25 Mq) that lead 
to the high probability region where the observed systems lie. 
Such high-mass donors in the LMBH models are relatively 
scarce. 

3.5. Effects of Radiation from the Accretion Disk 

Panel (d) of each figure shows the CMD of the donor stars 
alone. Comparing panels (c) and (d) in each figure, we see 
that the difference is most apparent in the IMBH scenario. 
Comparing panels (c) and (d) in Fig. 4, we see that the CMDs 
differ significantly only in regions corresponding to the giant- 
branch phases of the evolutionary tracks, and remain quite 



similar over the rest of the diagrams. This observation is sim- 
ilar for Fig. 5. And since the giant-branch phases lie in the low 
probabihty regions of the diagram, one can say that the optical 
flux from disk radiation contributes relatively little to the ob- 
served optical appearances of ULXs in the LMBH scenario. 
By contrast, for the IMBH model, a comparison of panels (c) 
and (d) in Fig. 6 shows that the contribution of disk irradiation 
and intrinsic disk radiation from viscous losses (see §2.5) is 
quite significant in this case, and has substantially altered the 
CMDs in most phases of the evolution. 

The greater influences of the disk radiation for the IMBH 
model (see Fig. 6) compared to those for the LMBH models 
seen in Figs. 4 and 5 can be explained by two effects. First, in 
general, the IMBH systems have higher luminosities by about 
an order of magnitude - largely by virtue of their higher mass 
donor stars. This obviously increases the power going into 
disk irradiation. However, a comparison of IMBH and LMBH 
systems at the same values of forb and Lpot reveals an impor- 
tant difference. The disks in IMBH systems still emit much 
more optical power than do the LMBH models, and this is 
predominantly at the blue end of the spectrum. The reason is 
that there is a substantial amount of intrinsic radiation from 
the disk (viscous release of gravitational potential energy; see 
§2.5). The physical explanation for this is straightforward. 
At the same radial distance from the central black hole, the 
IMBH system releases ~100 times more gravitational poten- 
tial energy than the corresponding LMBH system (simply due 
to the higher mass of the black hole). This release can end up 
dominating over the disk irradiation in regions where the bulk 
of the optical emission is released (see Fig. 2). 

Finally, it should be noted that the CMD for the donor star 
in a binary is expected to be different from that of a single star 
because of the mass loss via Roche-lobe overflow. 

3.6. Orbital Period versus X-Ray Luminosity 

Panel (b) in each of Figures 3, 4, 5, and 6 shows the evo- 
lution of Lpot with Porb- As can be seen in Fig. 3, the vari- 
ations in Porh reflect the changes in radius of the donor star 
through the various phases of evolution. Owing to the depen- 
dence of the Roche-lobe radius on the orbital separation, any 
rapid increase in the radius of the star while in contact with its 
Roche lobe causes a rapid increase in the orbital separation, 
and hence in Porb- The final rapid rise in Porb at the highest Lpot 
characterizes the expansion of the donor on the giant branch, 
leading to long orbital periods and high mass-transfer rates. 
As is apparent, the long values of Porb in the different mod- 
els indicate the giant-branch phases of the donors, and the 
short period regions pertain to the main-sequence and sub- 
giant phases of evolution. The intensity at a given location 
in panel (b) is proportional to the number of systems passing 
through the region and the amount of time spent by the sys- 
tems in that region. Consequently, the intensity in that region 
is a measure of the probability of finding a system at that lo- 
cation in parameter space. We then see that, for any model, 
the most probable regions for Lpot ^ 10^^ ergs s"' lie within a 
period range of ~ 1-10 days. 

It follows that the total probability of a system to be in the 
giant-branch phase is proportional to the sum of all the evo- 
lution times recorded in the long-period region, taken to be 
^orb ^ 10 days. Similarly, the total probability of being in the 
pre-giant phase is proportional to the sum of all the evolu- 
tion times recorded in the short-period region, i.e., Porb ^ 10 
days. Thus, we can calculate the relative probability for a 
ULX donor to be in the giant-branch phase versus that to be 
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Fig. 7. — The fraction of LMBH/IMBH binaries as a function of star cluster 
evolution time for three different lower limits on the X-ray luminosity. The 
different models are labeled in the panels. For each model, the blue curves 
represent fractions of the initial systems that have Lpot 10^^* ergs s"', the 
green curves represent fractions with Lpot ^ 2 X 10^^' ergs s"', i.e., ULX 
luminosities, and the red curves are for fractions of systems with Lpoi > iC" 



on the pre-giant branch by calculating the ratios of the two 
sums for all Lpot ^ 2 x lO''^ ergs s"'. The relative probabili- 
ties thus calculated are 0.08, 0.05 and 0.15 for Models La, Lb 
and Ic respectively. This reflects, in a quantitative way, the 
known fact that the systems are more likely to be found in the 
pre-giant phase than in the giant phase because of the rela- 
tively small amount of time they spend on the giant branch. 

The implication of the Poib-i'pot population diagram lies in 
the fact that given an observed value of Lpot, one can estimate 
the most probable orbital period corresponding to each of the 
three ULX scenarios presented in Figs. 4, 5 and 6. For in- 



stance, given an observed L, 



pot 



lO-*" ergs s-i. 



we see from 



Figs. 4, 5 and 6 that the most probable periods lie between 1- 
5 days, 1-5 days, and 2-10 days for the Models La, Lb and 
Ic respectively. Thus, the diagrams serve as rough guidelines 
showing what orbital periods are to be expected. 

3.7. Evolution of the X-Ray Luminosity 

Panel (a) of each of Figs. 4, 5 and 6 shows a population 
diagram of the evolution of the X-ray luminosity (Lpot) with 
the age of the system. Defining the ULX luminosities as 
^pot ?i 2 X lO-'^ ergs s"', we see that the systems in Model 
Ic are most active as ULX systems between ~2 and 15 Myr 
while systems in Models La and Lb are most active between 
~4 and 25 Myr These modest differences in active cluster 
lifetimes can be explained on the basis of the ranges of donor 
masses in the different models. Donors with higher masses 
have shorter nuclear evolution timescales, and hence shorter 
overall lifetimes. So, the higher the donor masses, the shorter 
will be the active lifetimes of the systems, leading to shorter 
overall timescales over which a star cluster containing such 
binaries is X-ray active. Shorter lifetimes of systems also 
mean that the mass-transfer rates are higher, leading to sys- 
tems with higher luminosities. Consequently, while systems 
in Model Ic can produce high Lpot systems even before as- 



cending the giant branch, systems with lower masses have to 
wait longer to produce their peak X-ray luminosities. It also 
follows that the numbers of systems having high Lpot will con- 
tain a higher proportion of higher-mass donors. 

Figure 7 elucidates this latter point. It shows plots of frac- 
tions of the initial systems that have Lpot greater than lO-'^, 
2 X lO^** and 10'*° ergs s"', at any given evolution time. We 
see that the maximum percentage of systems having (poten- 
tial) ULX luminosities at any time during the age of a star 
cluster is - 20% for Model La, - 30% for Model Lb, and 
~ 25% for Model Ic. We also see that the percentages of sys- 
tems having Lpot > 10"*° ergs s"^ are ~ 2%, 3%, and 20%, 
respectively, for these same models. It follows that the IMBH 
model allows for the largest numbers of very high luminos- 
ity systems. The percentages quoted above refer only to the 
most active phase of the corresponding clusters. The detailed 
evolution of the percentages with the star cluster age can be 
seen from the figure. However, the very small numbers of 
systems with Lpot > 10'**' ergs s"' at any time during the entire 
cluster lifetime make the LMBH models an unlikely scenario 
for the most luminous ULXs (independent of Eddington-limit 
arguments). 

3.8. Estimates of ULX Numbers 



Based on the L, 



pot - fev population diagrams, we can also 
make rough estimates of the numbers of ULXs of different 
Lpot that can be found in a galaxy in steady state, assuming 
the different models of ULX formation. 
Each element of the Lpot - fev matrix contains the number 

and 



of BH binary systems that have luminosities between Lpot ■ 
Lpot + ALpot and ages between fev and fev + Afev Both, Lpot 
and fev are logarithmically spaced with 100 bins per decade. 
Given the fact that we evolve 30,000 binary systems for each 
model, the fraction of all high-mass binary systems that lie 
in the above mentioned intervals of Lpot and fgy is given by 
m(Lpot,fev)/30,000, where m denotes the 700 x 700 Lpot -fev 
matrix. 

For the LMBH models, the BH accretors in this age in- 
terval, Afev, have all resulted from core-collapse supemovae 
(SNe) in an equal time interval at some epoch, fev, in the past. 
Assuming a uniform SN core-collapse rate of /?sn, the total 
number of BH binaries formed in this time interval is given 
by /?SN X /bh X Afev, where /bh is the fraction of all core- 
collapse supernovae that result in BH binaries with high mass 
donors (as determined from the BPS code). It then follows 
that the number of binary systems with luminosities between 
Lpot and Lpot + ALpot and evolutionary ages between fev and 
fev -I" Afev is given by: 

AA^(Lpot.j , fev.i) = X Afev,i X m(Lpot,j , fev,i) , (7) 

where j and / label the matrix element bins of Lpot and fev, 
respectively. The total number of sources with Lj > Lpot in 
a galaxy in steady-state can be obtained by summing over all 
the time bins of fev, and all the luminosity bins greater than 



-^pot? 



. s ^SN X /bh 

M> Lootk) = — X 

V pot,ky 30^000 

700 700 



(8) 



1 ^ev,i) 



/ ^ / ^ Afev.i X m(Lpot,j , iev,i 
j=k i=\ 

For illustration, we choose a uniform supernova rate 7?sn = 
0.01 yr"' (e.g., to represent a typical Milky- Way type galaxy). 
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Fig. 8. — The expected numbers of ULXs at the present epoch as a function 
of X-ray luminosity in a starburst galaxy with a core collapse supernova rate 
of 0.01 per year. The red, green, and blue curves correspond to Models La, 
Lb and Ic, respectively. 



and /bh ^ 3.75 x lO'^ and 9.44 x lO""* for models La and Lb, 
respectively. These values of /bh are deduced from the results 
of our BPS code. The computed steady-state cumulative lu- 
minosity functions for Models La and Lb are shown in the red 
and green curves of Fig. 8, respectively. 

For the IMBH model, we consider the IMBHs to be formed 
dynamically in young star clusters. Following a formulation 
similar to eq. (8) for LMBH models, we find for the IMBH 
model: 



^ysc X /iMBH X ./c 



cap 



30,000 

700 700 (9) 
^ ^ ^ ^ ^fev.i X m(Lpot,j , fev,i) j 

where, Ry^c is the formation rate of young star clusters (with 
^cius ^10* Mq) in a typical spiral galaxy, /imbh is the frac- 
tion of all such clusters that form an IMBH and /cap is the 
fraction of all such IMBHs that capture a massive companion 
into a close orbit. As illustrative values, we choose a uni- 
form /?ysc = 10"^ yr"' (estimated from the combined work of 
de Grijs et al. 2003; Hunter et al. 2003; Grimm et al. 2003), 
/imbh — 0.1, and /cap — 0.05 (see, e.g., Blecha et al.2006). 
The distribution of the numbers of systems as a function of 
Lpot is shown by the blue curve in Fig. 8. 

It can be seen from Fig. 8 that the expected numbers of ULX 
systems per galaxy are in fair agreement with the observa- 
tions (Ptak & Colbert 2004; Grimm et al. 2003) for the IMBH 
model, albeit with very large uncertainties. Specificially, 
these estimates in Fig. 8 indicate that a typical galaxy is ex- 
pected to harbor -0.07 and -0.03 ULXs with Lp^t > 2 x 10^^ 
and 10"*** ergs s"', respectively. By contrast, for the LMBH 
models, the expected numbers range from about —12-20 
for Lpot > 2 X lO"*^ ergs s"' per galaxy to —1.4-2.7 with 
Lpot > lO'*" ergs s"'. These numbers seem to be too large by 
about a factor of —100 compared with the observations (Ptak 
& Colbert 2004; Grimm et al. 2003). Moreover, there is con- 
siderably less flexibility, or uncertainty, in the LMBH models 
in comparison with the IMBH models. And, while one might 
dismiss the number of predicted ULXs with Lpot > 10"*° ergs 
s"' as violating the Eddington limit beyond the realm of the 
plausible, the predicted systems with Lpot > 2 x lO^*^ ergs s"' 
are not so easily dismissed. One completely straightforward 



way of reducing the computed numbers of ULX systems in 
the LMBH model is to adjust the A parameter downward (re- 
call that A is the dimensionless inverse binding energy of the 
stellar envelope of the BH progenitor). From the output of our 
BPS code, and an examination of Fig. 3 in Podsiadlowski et 
al. 2003, we can see that a value of A in the range 0.01-0.03 
will decrease the rate of production of BH binaries with high- 
mass donor stars by factors of — 10- 100. Such values of A are 
quite consistent with what stellar structure calculations tell us 
(see Fig. 1 of Podsiadlowski et al 2003; Dewi & Tauris 2000). 

Overall, there are a number of uncertainties associated with 
the various parameter values that we have used to calculate 
the ULX number estimates from eqns. (8) and (9). The pri- 
mary utility of Fig. 8 beyond providing crude estimates is that, 
given more accurate values for the uncertain input parameters, 
one can find the corresponding expected numbers of ULXs by 
simply scaling the expected numbers inferred from the figure. 

4. SUMMARY AND CONCLUSIONS 

In this paper, we have explored the optical properties and 
other observable parameters for a range of ULX models. We 
investigated a very large number of systems from three differ- 
ent ULX populations, two consisting of stellar-mass BH ac- 
cretors, which we refer to as stellar-mass (i.e., "low-mass") 
black holes (LMBHs), and one consisting of intermediate- 
mass black-hole accretors (IMBHs). For each population, we 
computed the evolution of 30,000 individual binary systems, 
and generated population diagrams to explore a wide region 
of parameter space. We computed optical CMDs for the donor 
stars and for the binary systems including radiation from the 
accretion disk. We also computed population diagrams show- 
ing the evolution of Lpot with donor age, and those showing 
the evolution of Lpot with Porb- 

When we plot the observed colors and magnitudes for six 
ULX systems found in the literature on the model CMDs, we 
find that all the data points lie in or near the high-probability 
regions spanned by the IMBH models on the CMD. On the 
other hand, fewer than half of the observed systems lie close 
to the high-probability region spanned by the LMBH mod- 
els La and Lb, respectively. In light of this observation, we 
conclude that IMBH models are somewhat favored by the op- 
tical observations. The locations of the data on the CMDs 
indicate that these systems correspond to high mass donors, 
with Mdon ^ 25 Mq. This provides an indication as to why 
the IMBH models are favored over LMBH models. Our BPS 
calculations for LMBH binaries do not produce very many 
successful systems with a donor mass of > 20 Mq, owing to 
limitations on the mass of the BH (—6- 15 Mq) and the re- 
quirement of stable mass transfer. On the other hand, the mass 
of the BH accretor in the IMBH models is set at 1000 Mq, al- 
lowing for very massive donor stars, and hence the ability to 
account for higher optical luminosities. 

We also find from the CMDs that the regions of high prob- 
ability correspond to the initial phases of evolution of the 
donors. This is apparent because all stars spend a predomi- 
nant fraction of their lifetimes on the main sequence. How- 
ever, this is particularly consistent with the color-magnitude 
observations, all of which lie in or very near the high proba- 
bility regions of the IMBH models. This leads us to conclude 
that the donor stars need to be predominantly on the main 
sequence or on the sub-giant branch. The B-V color range 
of about -0.35 to -0.10, corresponding to the high probabil- 
ity regions on the CMDs, indicates that the effective spectral 
class of the donors should be O through late B. We also dis- 
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cuss the effects of X-ray irradiation of the accretion disk, as 
well as intrinsic radiation generated viscously in the disk, on 
the different models. Disk radiation is found to be significant 
in the IMBH models, but less so in the LMBH models. We 
suggest that the main effects operating here are: (i) the higher 
donor masses in the IMBH systems, and hence higher values 
of Lpot, and (ii) the larger gravitational potential energy re- 
leased at a given radial distance in the accretion disk in the 
IMBH models. 

We compute population diagrams showing the evolution of 
Lpot with Poib- As concluded previously from the CMDs, the 
Porb-Lpot diagrams show that the giant branch phases of evo- 
lution in these systems form a low-probability region in pa- 
rameter space. It is seen that the most probable periods for 
all the models He in the range of ~1 - 10 days, corresponding 
to the main-sequence and sub-giant phases of evolution. We 
quantify the likelihood of ULXs being in the pre-giant phase 
by estimating the relative probabilities for the systems to be 
in a pre-giant phase versus the giant phase. We find the prob- 
abiUties to be in the ratios 0.08, 0.05 and 0.14 for Models La, 
Lb and Ic, respectively. We subsequently discuss the utiUty 
of the Porh-Lpot diagrams for estimating the probable range of 
forb for a system with a known Lpot. 

Finally, we compute population diagrams showing the evo- 
lution of Lpot with the evolution time of the host star cluster. 
Similar diagrams have been reported for various ULX models 
in our previous studies (see Rappaport et al. 2005; Madhusud- 
han et al. 2006). In the present paper, we use the diagrams to 
calculate the ULX production efficiency as a function of the 
age of a star cluster, for the different models. We also esti- 
mate the numbers of ULXs in steady state for a typical galaxy 
as a function of Lpot. We find, with crude estimates, that the 
IMBH models explain the observed numbers of ULXs rea- 
sonably well, with ~0.07 with Lpot > 10^^ ergs s"^ per Milky 
Way type galaxy to ^0.03 per galaxy with Lpot > 10'*" ergs 
s~^ For the LMBH models, the numbers range from about 
~12-20 for Lpot > 2 X lO^** ergs s"' per galaxy to ~1.4-2.7 
with Lpot > lO'*" ergs s~^ As discussed in §3.8 these num- 
bers are too large by about a factor of 100 compared with the 
observations (Ptak & Colbert 2004; Grimm et al. 2003). 

There are at least three ways of reducing the computed 
numbers of ULX systems in the LMBH model. One straight- 
forward way is to adjust the A parameter downward (see 
eq. [4] for a definition). We find that a value of A in the range 
0.01-0.03 will decrease the rate of production of BH binaries 
by the correct factors. Such small values of A are quite consis- 
tent with calculations of envelope binding energies of massive 
stars (see Fig. 1 of Podsiadlowski et al. 2003; Dewi & Tauris 
2000). Another possible explanation for the apparent overpro- 
duction of very luminous X-ray sources in the case of stellar- 
mass black-hole accretors is that the X-ray luminosities are, 
in fact, constrained to the Eddington limit, and that the large 
amounts of radiation-pressure ejected material severely atten- 
uate the soft X-rays coming from the vicinity of the black 
hole. Finally, we mention the possibility that for case B mass 
transfer (i.e., when the primary progenitors of the black holes 
are in the Hertzsprung gap during the time when mass transfer 
conmiences) the primary is more Hkely to produce a neutron 
star rather than a black hole even for fairly high initial masses 
(e.g., 50-60 Mq). As has first been argued by Brown et al. 



(1999) and then has been confirmed in detailed calculations 
by Brown et al. (2001), this has to do with differences dur- 
ing helium core buming. Basically the final structure of a 
star is very different depending on whether the star bums he- 
hum with a hydrogen-burning shell around it or without the 
hydrogen-burning shell (as one would expect if the star loses 
its envelope before or early during helium burning). With- 
out the hydrogen-buming shell the star ends up with a much 
smaller iron core at the end and most likely results in the pro- 
duction of a neutron star. This may be another reason why we 
do not find LMBH ULXs in nature. 

In estimating the numbers of ULXs, we have allowed for 
violation of the Eddington limit by using Lpot in our calcula- 
tions. This relaxation is imperative in order to explain ULXs 
using the LMBH models. However, the scenario by which 
the Eddington limit could be violated in a BH binary is not 
clear We note that the numbers for the LMBH models with 
Lpot > lO'*" ergs s"^ could be quite unphysical because of the 
fact that for these luminosities Lpot ^ LeiJci- The Eddington 
limit for a 20 Mq accretor is ^2.5 x 10^^ ergs s"'. Allowing 
for the fact that a massive donor star could be transferring he- 
lium onto the accretor during the later stages of its evolution, 
the Eddington limit is increased by a factor of 2, to ~5 x 10^^ 
ergs s"^ . Considering, however, that the bolometric luminos- 
ity may well be a factor of ~2-3 times the X-ray luminosities 
observed in the 1-10 kev X-ray band, the effective Edding- 
ton limit in the X-ray band is only ^1 .5 x lO-'^ ergs s"' . Thus, 
even if we allow for a violation of the Eddington hmit by rea- 
sonable factors of a few, and thereby achieve Lpot up to -^10''*' 
ergs s"', it is difficult to understand X-ray-band luminosities 
above ^5 x lO-'^ ergs s~'. 

As for the IMBH model, we have not studied in any de- 
tail the formation scenarios of IMBHs or of the capture of 
binary companions by the IMBHs (but see Portegies Zwart 
et al. 2004; Tutukov & Fedorova 2005; Blecha et al. 2006). 
Despite much evidence in the literature supporting IMBH bi- 
naries as candidates for ULXs, the formation mechanisms of 
IMBHs are not clear. In particular, it is unclear how supermas- 
sive stars evolve, and whether they can undergo core collapse 
to form an IMBH. In calculating the numbers of ULXs in this 
scenario, we have assumed plausible values for the formation 
rate of young star clusters in a galaxy, the fraction of these 
clusters that successfully produce IMBHs, and the probabihty 
of capture of a massive companion by the IMBH. 

Finally, we should keep in mind the possibility that neither 
the IMBH nor LMBH binary models may be correct in ex- 
plaining the most luminous ULXs. It is conceivable that the 
most luminous ULXs may, in fact, be free floating IMBHs that 
are fed by tidally disrupted passing field stars (see Volonteri 
& Pema 2005, for a related scenario). 
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